clear;clc;

load csp.mat;
load cspd.mat;

cspQ=csp-cspd;

load acoustic\csp.mat;
load acoustic\cspd.mat;

cspA=csp-cspd;

cspAQ=cspA;
cspAQ(:,401:800)=cspQ(:,401:800);

[zz,xx]=size(cspAQ);
dx=5;
dt=0.0005;
z_tick=[0:1000:zz];z_label=dt*z_tick;
x_tick=[0:100:xx];x_label=dx*x_tick;
Font=14;

plotimage(cspAQ);
set(gcf,'units','normalized','position',[-0.1 0 0.5 2.0]);
ylim([0,zz]);xlim([0,xx]);
set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x(m)','FontSize',Font);
ylabel('Time(s)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);
